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^ ; 1 Introduction 

CLc In this paper we study form factors of quantum integrable models in finite volume. The possibility 
to obtain exact non-pertubative results in this field is of great importance. The Quantum Inverse 
Scattering Method created by L. D. Faddeev and his school fl], 0, f| allows one to formulate general 
methods of solving this problem. Namely, in order to obtain an explicit representation for a form 
factor of an operator O, it is sufficient to embed this operator into the algebra 

^; R(\,fi)[T(\)®T(ji)) = \T(ji)®T(\))R(\,n). (1.1) 

After this is done, the problem is reduced to the standard Algebraic Bethe Ansatz procedures. 

We would like to mention that our approach is drastically different from the ones, used in 
relativistic models ||]-[[n]]. We consider form factors between the states constructed on the bare 
vacuum, instead of the physical one. We also do not use the axioms, describing analytical prop- 



erties of form factors. All our calculations are based on the algebra (Id) only. In spite of this 
approach meets serious combinatorial difficulties, nevertheless it provides us with explicit deter- 
minant representations for form factors. 

In the present paper we consider Quantum Nonlinear Schrodinger equation (QNLS) and XX Z 
(XXX) Heisenberg chains of spin 1/2. The QNLS model in a finite volume L with periodically 
boundary conditions is described by the Hamiltonian 

#qnls = dx (d x tfd x V + , (1.2) 
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where ty(x,t) and W(x,t) are canonical Bose fields 

[*(x,t),tf(y,t)]=6(x-y), *|0) = 0, (0|^ = 0, (1.3) 

and c is the coupling constant. 

The Hamiltonian of the XX Z Heisenberg chain with M sites is 

M 

H xxz = -J2 (^- m) ^l m+1) + 4 m) ^- m+1) + A(«r^VJ m+1 J - 1)) . (1.4) 



m=l 



Here cx^ m ^ are local spin operators (Pauli matrices), associated to the m-th site, a± = \{o x ± % 



i a 



yji 



a (M+i) _ ^.(i)^ Parameter A defines the anisotropy of the model. Particular choice A = 1 
corresponds to the XXX chain. 

Let us describe now the content of the paper. In the section 2 we recall the basic definitions of 
the Algebraic Bethe Ansatz. In the section 3 the determinant representation for the form factor 
of the local field \I/ of QNLS is obtained. In the section 4 we consider different representations for 
form factors of QNLS and prove their equivalency. The last section is devoted to the form factors 
of Heisenberg chains. 

2 Algebraic Bethe ansatz 

We consider integrable models, which can be solved via the Algebraic Bethe Ansatz. Let the 
monodromy matrix is 

and the highest vector |0) (pseudovacuum) of this matrix has the following properties: 

A(A)|0> = a(A)|0>, D(A)|0) = d(A)|0), C(A)|0> = 0. (2.2) 

The dual vector (0| is defined similarly 

(0L4(A) = a(A)<0|, (0|D(A) = d(A)<0|, (0\B(X) = 0. (2.3) 



The explicit form of functions a(A) and d(X), entering Eqs. ( |2.2| ), ( |2.3| ), depends on specific model. 

The i?-matrix, defining the commutation relations of the monodromy matrix entries ( |1 . 1[ ) , can 
be written in the form: 



( f(fjL,X) \ 

<?(//, A) 1 

1 g(ji,X) 

V f(ji,X)J 



(2.4) 
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For the models QNLS and XXX Heisenberg chain the i?-matrix is rational: 

/(ftA) = /^A+i£ 9 (,,A) = ^. (2.5) 

Here c is a constant. Below we shall use also two new functions 

w >n _ fifrX) fi-X + ic 

h{V,X) =—, — rr = : , 

g{H,\) ic 

(2.6) 

t(n A) = 9 ^ ,X ^ = 

/i(/z, A) (/i — A)(/i — A + ic) 

For the XX Z chain the i?-matrix has the same form (|2.4|) , however the functions g(A,/i) and 
f(\,n) become trigonometric ones 

sinhf/i - A + 27) , ishi7 

smh(/i — A) smh(/i — A) 

Respectively the functions h(X, //) and t(X, fi) should be replaced by their trigonometric analogs 
according to the definition (|2.6|) . The constant 7 is related to the anisotropy of the model by 
cos 7 = A. 

The entries of the monodromy matrix ( pTT ) act in a space, consisting of states 

\v N ({X})) = l[B(x j )\o), N = 0,1,... (2.8) 

Here {A} are arbitrary complex parameters. Similarly dual states can be constructed by the 
operators C(A) 

N 

PM{A})| = <0|nC(Ai), AT = 0,1,... (2.9) 

i=i 

In what follows we shall refer to the states ( |2.8|) , ( P-9|) with arbitrary A-s as £ BA- vectors' (Bethe 
Ansatz vectors). 

The transfermatrix r(A) = tr T(A) = A(X)+D(X) generates the complete set of the conservation 
laws of the model. The eigenstates of the transfermatrix have the form ( |2.8| ), ( |2.9| ), however the 
parameters {A} are not arbitrary, but satisfy the system of Bethe equations 

•™n$£g = l. where r(A) = $J. (2,0, 

The eigenvalues of the transfermatrix are 

rO*)IM{A})) = e(A*|{A})|^({A})), 

TV jv (2.11) 

0(/i|{A}) = a(//) J] f{v,Xm) + d{v) II 

m=l m=l 
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Let us focus now at the models, which will be considered in the next sections. For the QNLS 
the pseudovacuum (|2.2|) coincides with the Fock vacuum (|1.3|) . The coupling constant c is equal 
to the constant in ( [2.5|) , the functions a(A) and d(X) are 



a(A) = exp(-iAL/2), d(X) = exp(zAL/2). (2.12) 

In the XX Z chain with M sites the pseudovacuum vectors |0) and (0| are 

|0> = ® Tm, (0| = (|0)) T . (2.13) 

m=l 

For intermediate calculations the inhomogeneous XX Z chain is suitable. For this case the mon- 
odromy matrix is given by 

T(A) = L M (X - £ M ) • • -Li(A - 6), (2-14) 
where £ m are arbitrary complex numbers, and local L-operators are 



cosh ( A — at m ) s i n ^<j[ m ) 



L m {X)=t\ \ fa' > , ' s (2.15) 

sm74 ' cosh( A + z^er™ J / 

The functions a(A) and d(X) are equal to 

M M 

a(A) = J] cosh(A - ^ m - 17/2), d(X) = J[ cosh(A - £ m + ij/2). (2.16) 

m=l m=l 

In the limit £ TO — > we turn back to the homogeneous XXZ chain. 

Particular case A = 1 corresponds to XXX Heisenberg magnet. This model is described by 
rational i?-matrix ( [2.4|) with c = 1. The vacuum eigenvalues are 



M M 

a(A) = J] ^ - ^ - i/2), d(A) = I] ( A - ^ + */2)- (2.17) 

m=l m=l 

In the conclusion of this section we would like to mention that the Hamiltonians, presented in 
the Introduction, are not of the most general form. One can consider, for instance, Heisenberg 
magnets in external magnetic field. Similarly it is possible to consider QNLS with a chemical 
potential. In order to do this, one need to add to the Hamiltonians ( |1.2| ), ( |1.4| ) the corresponding 
terms. However these additional terms do not play an essential role in the models in finite volume. 
They become significant only in the thermodynamic limit, in particular they define the structure 
of the ground state and the thermodynamic properties of the models. Since the thermodynamic 
limit is not considered in this paper, therefore we restrict our selves with the Hamiltonians {T72) 
and (|TJD only. 
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3 Form factor of the local field in the QNLS 



The calculation of form factors in the framework of the Algebraic Bethe Ansatz formally can be 
reduced to the calculation of scalar products of a special type: 

s N m,w) = (oi m <?(^)j (n^)j io>. (3.i) 

In the last equation at least one of the states, say (0| IL=i C(fij), is an eigenfunction of the 
transfermatrix. In other words the parameters fix, . . . , //jv satisfy the system of the Bethe equations 
(|2.10|) , while the parameters Ai, . . . , Ajy remain free, i.e. Ylf=xB(X j ) is an arbitrary BA-vector. 

Let us demonstrate the general idea of the method, which will be used. Consider a form factor 
of an operator O 

' N v 




n^,)]|0). (3.2) 

Here both states (0\lff=i CQij) and Y\f=i B (vj)\0) are ekrenstates of the transfermatrix. The 
operator O acting on the vector B(uj)\0) gives a new state, which generally speaking is not 
an eigenstate of t(A). However, if this new state can be presented as a linear combination of 
BA- vectors with numerical coefficients 





|0> = L'V. inWMlO). (3.3) 
then evidently the form factor Fq is equal 

*b = £c fc <0| (nC(^)) (n^CAf)] |0>. (3.4) 
k \j=l J \j=l J 

Anyway the state (0| ITj=i C(a*j) remains the eigenstate of the transfermatrix. Thus, in order to 
find explicit expression for the form factor O, one need to calculate the scalar products (|3.1| ). In 
order to do this, it is enough to use the algebra ( p. . 1|) only. Moreover the specific form of the vacuum 
eigenvalues a(A) and d(X), as well as the structure of the pseudovacuum, are not essential. However, 
the straightforward calculation of the scalar products appears to be a complicated combinatorial 
problem. These difficulties had been partly overcame in ||16|| , [ jT7|| . In the paper |18j the determinant 



representation for the scalar product of two arbitrary BA- vectors was obtained in terms of auxiliary 
quantum operators — dual fields. The explicit determinant representation for the scalar products 
(O) was found in m (see also EDI): 



N °\ I N 



s N m, {A}) = n d(»a)d(\ a ) n k^, a 6 



x < II 9{Xa, h)g{fib, Ha) \ det (M jh ) , (3.5) 

N>a>b>l 
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where 

M jk = t(ji h , A,) - r(A J )t(A J , /*) ft t^tt - (3.6) 

Recall that here the parameters {//} satisfy the Bethe equations ( |2.1U| ), while the parameters {A} 
are free. Before applying this formula to the calculations of form factors we would like to mention 
several features of the representation ( |3.5|) . 

This representation is valid for arbitrary model, possessing rational i?-matrix ( |2.4|) or trigono- 
metric one. The representation ( |3.5| ) also allows one to prove the orthogonality of the transferma- 
trix eigenstates directly. Indeed, the particular case of an BA- vector nJLi B(Xj)\Q) is an eigenstate 
of r(A). Therefore, one can demand the parameters {A} to satisfy the Bethe equations ( |2.10 ), but 
to be different from the parameters {/i}. In this case we obtain the scalar product of two different 
eigenstates, which must be equal to zero. It is easy to check, that for the such choice of {A} the 
determinant detM in ( |3.5| ) does vanishes. Indeed, the matrix (STB) turns into 

M jk = t(n k , Xj) + Vjt(Xj, ii k ) (3.7) 

where 

V- = TT ^m)h(X m , Xj) ^ ^ 

m=1 h(ji m , Xj)h(Xj, X m ) 
This matrix has the eigenvector £ with zero eigenvalue: 

N 

n g(jJ-k,Pm) 

N m=l 

J2 MjA = 0, = ^ . (3.9) 

k=1 n g{^k,X m ) 

m=l 

The proof of the ( |3.9| ) is given in Appendix 1. Hence det M = 0. 

Now let us apply the representation (|3.5| ) to the calculation of specific form factors. Namely, 
consider the form factor of the local field ty(x,t) in the QNLS 

F 9 {x,t) = (0| ( nC(^)j f n^^oj l«>- ( 3 - 10 ) 

The dependency on time t and distance x evidently separated 

!N N+l "I 

(itfil — ix/i/A — ^ (^A| — ixXj^j > F^(0, 0). (3.11) 
k=l ' j=l I 



The action of the operator \l/(0, 0) on the BA- vector had been found in ITS 



N+l N+l / N+l \ N+l 

*(0,0) n^)|0) = -2v^E a (^) Uf^t^m)] IT ^(A m )|0>. (3.12) 



3=1 1=1 



m= 1 / tn= 1 
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In spite of the parameters Ai, . . . , Xn+i satisfy the Bethe equations 



r(A,) = I] 1/l(Ap ' Aj 



(3.13) 



p =i h(\j, A p ) ' 

the states entering the r.h.s. of (3.12) are not eigenstates of the transfermatrix, since each of them 
depends only on N parameters A. Nevertheless all these states are BA-vectors, therefore we can 
apply fl3.5l ) for calculation of the scalar products 



JV 



AT+l 



(3.14) 



since the state (0| Uj=i Ojij) remains an eigenstate. After simple algebra BTJ we obtain for the 
form factor F\y(0, 0): 

N N+l N+l 

F«(0,0) = -iy/cY[g(/j, a ,ij, b ) Y[g(X b ,\a) Yl h (^a, X b ) 

a>b a>b a, 6=1 



X 



N N+l /N+l 

nwn^) 



k 



a=l 



6=1 



Here we introduce new matrix S: 



N 

n K^Xj) 

Sjk = t(fik, Aj)^— j t(Xj,fi k ) 



\£=1 



N 



N+l 

P =i 



n h(Xj,ii m ) 

m=l 
N+l 

n h(Xj,x p ) 
P =i 



3,k = 1, 



(3.15) 



AT. 



(3.16) 



The matrix 5^ can be obtained from the matrix S via the replacement of the £-th row by S 



Sjk, 



The sum with respect to £ in (|3.16|) gives a single determinant: 



AT+l 

^ (-l) <+1 detS$ = det(5- ifc - S N+1>k ). 



N+l,k- 



(3.17) 



(3.18) 



In order to prove (|3.18|) one can use the Laplace formula for determinant of sum of two matrices, 
and to take into account that S'jv+i.it is the first-rank matrix. Thus we arrive at the determinant 
representation for the form factor of the local field ib QNLS ||21|| : 



N N+l N+l 

Fvp(0,0) = -i^Y[g(fi a ,fi b ) Y[g(X b ,X a ) J[ h(X a ,X b ) 

a>b a>b a, 6=1 



N N+l 

x ]J II d(X b )det(S jk - S N+1>k ). 



(3.19) 



a=l 



6=1 
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Here Sjk is given by fl3.16|) , and Sjv+i,fe can be obtained from Sjk via the replacement Xj — > Ajv+i- 
In the conclusion of this section we would like to emphasize once more that the most serious 
problem of the described method is to embed the operator O, whose form factor wanted to be found, 
into the algebra ( TTT|) . However, after the commutation relations between O and the monodromy 
matrix are defined, then the problem can be solved via standard methods of the Algebraic Bethe 
Ansatz. 



4 Other representations for form factors 

As we have mentioned already, the method used for calculations of form factors of relativistic 
models is essentially different from the method, presented in the previous section. Nevertheless 
a series of the results, obtained in this domain for the quantum field theory models, also can 
be applied for the models in finite volume. Thus, for instant, the form factors of 0(3)-invariant 
cr-model unsuspectingly had been found closely related to the form factors of the QNLS p^| . In 
the present section we consider several examples. 

The approach, developed in papers |0| [XT], allows to express form factors of the QNLS in 
terms of some rational functions S a (originally they were denoted as a a ), which depend on iV M 
arguments {/i a } and N\ arguments {Aft}. The only singularities of these functions are simple poles 
at Xj = fik, and the residues in these points can be expressed in terms of S a , depending on the 
arguments {/i a -^} and {A^}. This property allows one to construct S a via a recurrence, however 
during a long time this recurrence have not been solved in general form. 

On the other hand the same functions S° appeared to be useful for the analysis of the form 
factors of 0(3)-invariant a-model. The series of the explicit determinant representations for E a , 
depending on elementary symmetric polynomials er| n , were found in |23 

1 Q{n-k) 



a 



(n) 



{x\ , . . . , x n ) 



(n — k)\ &r( n " fc ) w , 
The explicit solution for the functions S a has the form: 

s a (M,{A}) <l0tO 



Y[ (x + x m ) 



2=0 



A(ji) A(A) rL 6 0"a - A 6 )' 



(4.1) 



(4.2) 



Here (N„ + A^^) x (N^ + N\) matrix Q has the following entries: 



M, 



a 



(Nx+N^-l) 



{fia^k - —}, {X a + — } 



-e a af x+N * 



IC, 



{^a^k + — }, {A a 



2C- 



N, 



(4.3) 



A; 



ok 



a 



\Va + — },{A a ^fc - — }J 
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—a 



(N x +N^l) 
3 



{^a ~ y}, {K^k + , fc = l, 



(4.4) 



Index j runs through 0, 1, ... , iV M + N\ — 1. By A(-) we denote the Wandermonde determinants. 

Form factors of the QNLS can be expressed in terms of the functions S a . In particular the 
form factor of the local field F^(0,0), considered in the previous section, is proportional to the 

A 



function S° with = iV A - 1 = N: 



N N+l 

F*(0, 0) = -zVB J] d(fxa) II d ( X b) S a (Miv, Wn+i] 

a=l 6=1 



a=0 



(4.5) 



It is wort mentioning, that in this case Q is (2N + 1) x (2iV + 1) matrix. On the other hand, as we 
have seen above, the form factor of the local field is proportional to the determinant of the N x N 



matrix fl3.6|) . There exists a nice method to reduce the determinant of the matrix Q (|4.3| ), 
to the determinant of (|3.6|) . We shall demonstrate this method for more simple form factor. 

Consider an operator Qi of the number of particles in the interval [0,x\. In the QNLS this 
operator is equal to 

(4.6) 



Q x = r¥(y)*(y)dy. 
Jo 



The form factor of this operator was calculated in |l7j in terms of the function S a with 

N X = N: 



N 



N 



F Ql = {0\[l[C( f , j )]Q 1 [l[B(X j )\\0) 



exp i.ix Y^i^a - A*o)| - 1 



^({»} N ,{\} N ) 



(4.7) 



a=0 



Let us show, how one can reduce the determinant of 2N x 2N matrix ( |4.2|) to a determinant of 
N x N matrix. 

Theorem 4.1 The function E a ({/x}jv, {A}n) is proportional to the determinant of N x N matrix 

N N 

£°({/i}v,{A}Ar) = ng(\j,\ k )g(vk,Vj) II h (^a,Vb) 

a>b a, 6=1 



x det (e a Vjt(fi k , Xj) + t(X j} n k )) , (4. 

where 

N 

v.= T t t . 

=1 ^(Aj, Hm)h(X m , Xj) 



y _ -Q %«; Aj)/l(Aj, A m ) ^ ^ 

m= 



(compare with 
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Proof. The determinant representation for the function T, a ({fi}^, {X}n) is given by 
with iV M = N\ = N. Let us multiply the matrix Q from the left by matrix U 

U hk = P jk = vf~ k '\ J = 1,...,N, (4.10) 



CW = Qok = qf' k -\ j = i,...,N. (4.11) 

Index k runs through the set 0, 1, ... , 2N — 1, parameters pj and qj are some complex numbers, 
which will be fixed later. Obviously 

and 

TV 

det C/ = A(p)A( 9 ) J] (P« - ft)- (4-13) 

o,6=l 

The product f/fi can be written as 2 x 2 block-matrix 

where each of these blocks is x matrix. Using the representation for elementary symmetric 
polynomials ( |4.1| ), one can easily find the entries of these blocks, for example: 

(PM), fe (M, m, {A}) = n (pj + a*. - 1) n (ps + A « + f ) 

a^k 

-e a n (pi + a*. + 1) n (pi + A « - f ) > ( 4 - is ) 

Other blocks have similar form. Namely, replacing pj by qj in ( |4.15| ) we obtain the block QM. 
As for the right blocks, they can be obtained from the left ones, in which one should replace 
{fi} by {A} and put a = 0. Thus we obtain the representation for the function T, a ({fi} N , {X}n), 
containing 2N arbitrary complex parameters pj and qj, but does not depending on their specific 
choice. 

Using well known formula for a determinant of a block-matrix 

det ( C D ) = det A ' det ^ D ~ CA ~ lB )> ( 416 ) 

one can reduce det(UQ) to the determinants of N x N matrices. However, in order to simplify 
the calculations, we shall fix the parameters pj and qj. Let, for example, pj = Xj — ic/2 and 
qj = Xj + ic/2. Then we have for two left blocks: 

N N 

(PM) jk = -e a l[^ a - X 3 ) U(X aj - ic), (4.17) 

a =! a=l 

a^k 
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N N 



(QM) Jk = no*«- A i)ii( A «w+ ic )> ( 4 - is ) 

Here \ a j = A a — Xj. Observe that 



a=i a=l 



(QM) jk = - Uj (PM) jk , where Uj = e~ a ]\ ( Xa *+ lc \ , (4.19) 

0=1 lc ) 

Hence, due to ( [OB ) we have 

det(Utt) = det(PM) jk ■ det((QA) jfc + u j (PA) jk ). (4.20) 
After simple algebra we arrive at 

det(Un) = WdetG (4.21) 



where 



and 



JV N 

W=U [(// a - X b ){\ ab -ic){\ ab + ic)} n^aftAto], (4.22) 

a,b=l a>b 



N 

n A aj 

^ = ^-^3 — je« n ( fi \- Xj i +ic ) - n ilLa ~ Aj ~ ^ 



v 

II 0*a - Aj) 



a=l 



a =i V A a j + zc y a=1 y A a j zc 



In order to reduce detG to the determinant (|4.9|), one need to make one more step: 



where 

N 

1 n(Aj-/i a ) A 

I> = Q=1 N , and detr=n(^- (4-25) 

j ^ k n a a a>b 



The product GT can be computed explicitly via the method, describing in the Appendix 1. The 
resulting representation for the form factor Fq 1 exactly coincides with (|4.8| ) (|4.9| ). 

One can use the same idea in order to reduce the determinant of (2 A" + 1) x (2 A" + 1) matrix 



to the determinant of the N x N matrix (3J3). A small complication is caused by the fact 



that in this case the blocks in the matrix ( |4.14j ) have different dimension, for example, the block 
PM is N x N matrix, while the block QM is (iV + 1) x A" matrix etc. This leads to the appearance 
of the term Sjv+i,& in ( |3.19| ). 
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Thus we have demonstrated different representations for form factors and proved their equiva- 
lency. The representations in terms of elementary symmetric polynomials are quite characteristic 
for the form factors of the relativistic models (see, for instance, ||10|| ). At the same time the rep- 



resentations (|3.19|) , ( |4.8|) look more natural, since they depend only on the functions, entering 
the i?-matrix of the model. These representations also are suitable for the computation of the 
correlation functions. 



5 Form factors of Heisenberg chains 

As we have seen already, the calculation of the form factors is based on the explicit determinant 
representation (^.5[), ( |3.6| ) for the scalar products of the special type (|3.1p. However the calcu- 
lation it self of the scalar products in the framework of the Algebraic Bethe Ansatz deals with 
rather complicated combinatorial problems. The method based on the explicit representations 
for the functions S a in terms of elementary symmetric polynomials also makes use of solving of 



complicated recurrence. In the paper [24] a new approach was proposed, which allows to simplify 
significantly the calculation of scalar products and form factors. There a factorizing Drinfel'd twist 
F was constructed. The new basis generated by the twist F is extremely suitable to study the 
structure of BA-vectors. and for representation of local operators in terms of the entries of the 



monodromy matrix. In the paper [23 this method was used for calculation of the form factors in 
Heisenberg chains. We would like to present here the main results of the paper [EH] . 
Consider inhomogeneous XX Z (XXX) model ( [2.14|) : 



T{\) = L M {\-£ M )---L 1 {\-Z 1 ). (5.1) 
Then the local spin operators can be presented in terms of the monodromy matrix entries as 



Am) 



/m-l \ I M 

n^.)p(*m) n rfco 

\a=l / \a=m+l 



/m-l \ / M 



0+ ' = I II )C(U[ LI (5.2) 

\a=l / \a=m+l 

/m-l \ I M 

° { z m) = n^a) )(Mu-D(u)[ n <^ 

\a=l / \a=m+l 

Here r(A) = A(X) + D(\) is the transfermatrix. 

The representation ( |5.2| ) allows to find the form factors of local spin operators of XX Z (XXX) 
chains. Consider, for example, form factor F^ of the operator a^: 

(N+l \ ( N \ 

fH* = (0| I n C(nM f J] B(\ 3 ) \ |0). (5.3) 
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Both states in ( |5.3j ) are eigenstates of the transfermatrix, hence substituting here ( |5.2| ) we obtain 

m-l M /N+1 \ / TV \ 

f m = n ^im) n w) - <oi n ^ ^(u n io>. (5.4) 

a=l o=m+l \j=l / \i=l / 



Here 9 is eigenvalue of the transfermatrix ( 2.11 ). Thus we again have reduced the form factor to the 
scalar product of eigenstate by BA- vector. It is enough now to use the representation fl3.5|) , which 
gives us the representation for the form factor F ( m ) in terms of a determinant of (iV + 1) x (N + 1) 
matrix. 

Let us present the final result in the form given in [p5|l . Introduce (N + 1) x (N + 1) matrix 

d 

Tjk({v}N+i\{v}N+i) = —9{^\{v} N+1 ). (5.5) 

Ol>k 

Then the form factor is proportional to the ratio of two determinants 



(m) _ Ua=i Ub^i fiva, 6) detr^MK^A!,...^^}) 



(5.6) 



nli nr=i /(Ao, 6) det^Mi^A!,...^^})' 

Here i^-fc is the Cauchy matrix 

K 3k({n}\{v}) = V^i^i - ^fe), (5-7) 

where <p{X) = A for XXX chain and y(A) = sinhA for XX Z chain. Actually the function y) _1 (A) 
coincides with the function g(X) up to constant factor. The representation similar to ( |5.6|) also 
exists for the scalar products (K 



6 Conclusion 

We considered form factors of exactly solvable (completely integrable) models in the finite volume. 
We demonstrated that the form factors can be represented as determinants of matrices. The 
dimension of the matrix is related to the number of particles in the corresponding state. In 
our earlier publications we showed that this representation is useful for the theory of correlation 
functions. The contributions of all form factors can be taken into account and the correlation 
function also can be represented clS db determinant. 



A Zero eigenvector 

Consider for simplicity the case of the rational i?-matrix. Let N x N matrix ( |3.7|) is given 

M jk = t{fi k , Xj) + Vjt(\j, fi k ). (A.l) 
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Here 

V = f] ( A i - ^ m + ic)(X m - Xj+ic) 2 ^ 
3 mil (Vm ~ Xj + ic)(Xj - X m +ic)' 
Recall also that t(X, //) = (ic) 2 / (A — /i)(A — /i + ic). Let us prove that 

N 

n n (i^k — A m ) 

X)M ifc & = 0, where & = ^ ■ ( A ' 3 ) 



fe=i 



m — 1 



We have 



where 



iV 

£ MiifcCfc = («c) 2 ( G + + V$G_ ) , (A.4) 



k=l 



N 

n y n {y-k — x m ) 

G± ~ S (^-A.O^-A.-iic) ' g (/ifc _ /im) ' (A ' 5) 



m— 1 



In order to find G± consider auxiliary integral 

j ± = J_ / d A T7 ~ Am ) (A 6) 

± 2m 7 (z - A 3 -)(z - Aj ± ic) ^ (z - fi m ) ' 

The integral is taken with respect to the closed contour around infinity (in other words all the 
poles of the integrand lie within the contour). Obviously the integral I± is equal to the residue at 
infinity, which in turn is equal to zero, hence I± = 0. On the other hand the value of the integral 
I± is equal to the sum of the residues in the poles z = fi k , k = 1, . . . , N and z = Xj — ic. Thus, 
we obtain 

4 = C ±T in(^). (A.7) 

IC£=1 \Vm-X J ±lcJ 

The first term in the r.h.s. of (|A.7|) corresponds to the residues at the points z = fik, k — 1, . . . , N; 
the second term corresponds to the residue at the point z = Xj — ic. Hence 

G± = ±1 ft f A "-^f' C ) . (A.8) 



Substituting G± into ([A.4 ) and using explicit expression (|A.2|) for Vj, we arrive at 



N 



J2 M jk £ k = 0. (A.9) 



fc=i 



In the case of trigonometric i?-matrix, after the replacement exp{2A} — > X and exp{2/i} — > // 
the proof reduces to the rational case. 
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